submitted to PRA 



A versatile electrostatic trap 



O 
O 



Jacqueline van Veldhoven^'^, Hcndrick L. Bcthlcm^''^, Melanie Schnell^, and Gerard Meijer^ 
^ Fritz-Haber- Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany 
^ FOM-Institute for Plasmaphysics Rijnhuizen, P.O. Box 1207, NL-3430 BE Nieuwegein, The Netherlands 
Laser Centre Vrije Universiteit Amsterdam, de Boelelaan 1081, NL-1081 HV Amsterdam, The Netherlands 

(Dated: February 2, 2008) 

A four electrode electrostatic trap geometry is demonstrated that can be used to combine a dipóle, 
quadrupole and hexapole field. A cold packet of ^^NDs molecules is confined in both a purely 
quadrupolar and hexapolar trapping field and additionally, a dipóle field is added to a hexapole 
field to create either a double-well or a donut-shaped trapping field. The profile of the ^^NDs packet 
in each of these four trapping potentials is measured, and the dependence of the well-separation and 
barrier height of the double-well and donut potential on the hexapole and dipóle term are discussed. 



PACS numbers: 33.80.Ps, 33.55.Be, 39.10.+j 
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I. INTRODUCTION 

The creation of cold molecules offers unique possibili- 
ties for studyinff cold coUisions dipole-dipole 
interaction ^ Ei > doing high-resolution spec- 

troscopy H, 13 . In each of these fields, traps are either 
necessary, or they enhance the potential of the method 
considerably. Electrostatic traps in particular are highly 
versatile. They are rather deep (on the order of 1 K) and 
by changing electrode geometries, many different trap- 
ping potentials can be formed Efl El lil El El El- 
Here we demonstrate an electrode geometry that is used 
to forni different trapping fields by changing the voltage 
applied to the electrodes. In the trap a predominantly 
dipolar, quadrupolar or hexapolar field can be created. 
When implementing the latter two, a pure quadrupole or 
a pure hexapole trapping potential is generated. Combin- 
ing a dipóle field with a hexapole field gives rise to two 
additional trapping potentials, a double-well or a ring- 
shaped potential, depending on the sign of the dipóle 
term. As voltages can be switched on relatively short 
timescales, rapid changing between these trapping fields 
is feasible as well. 

A double-well trap in particular has proven to be fruit- 
ful in atomic physics experiments. For instance, double- 
well potentials have been used in interference experi- 
ments to study coherence properties of BECs E3 and to 
visualize phase singularities associated with vortices in a 
Na BEC 17]. Furthermore, two weakly coupled BECs 
in a double-well potential have been shown to form a 
Josephson junction and coUision studies have been 
performed in a double-well interferometer, enabling the 
detcrmination of the s- and d-wave scattering amplitudes 
of «^Rb 19]. 

In electrostatic traps the formation of a double-well po- 
tential with variable barrier height and well separation is 
relatively straightforward. Along with the possibility of 
rapid switching to other types of trapping fields, this cre- 
ates favorable conditions for the study of cold coUisions. 
Molecules trapped in either of the two wells of a double- 
well trap will gain potential energy when switching to 
a single- well trap, resulting in two packets of molecules 
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FIG. 1: Scheme of the trap. The trap is cylindrically sym- 
metric and consists of two parallel rings in the middle and 
two endcaps at the sides. 



that are accelerated towards each other. As they reach 
the center of the trap, they will coUide with a velocity 
that is dependent on the shape of both the single- and 
double-well potential. As these shapes can be readily 
varied, the velocity of the molecules is tunable. Combin- 
ing this feature with an inherently low velocity spread 
of the trapped molecules makes rapid single- to double- 
well conversion an excellent starting point for studying 
coUisions at low temperatures. 



II. THE VARIOUS TRAPPING 
CONFIGURATIONS 



In a cylindrically symmetric gen rnetr 
tential <Í>(p, z) can be expressed as [lil i 
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with zq a characteristic length scale. Here, the first 
term is a constant voltage, the second term a constant 
electric field, the third a quadrupolar electric field and 
the fourth a hexapolar electric field. A molecule with 
a quadratic Stark shift will experience a harmonie force 
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in a quadrupolar electric field, whereas a hexapolar field 
is needed to form a harmonie trap for a molecule with 
a linear Stark shift. In our experiments we use ^^NDs 
molecules in the low-field seeking hyperfine levels of the 
\J,K) — |1,1) inversion doublet of the vibrational and 
electronic ground statě. These molecules experience a 
quadratic Stark shift for small electric field strengths and 
a linear Stark shift for larger electric field strengths. The 
dependence of the Stark shift on the electric field strength 
E for ^^NDs in these levels is given by pll |: 

Ws.ar.{E) = + - , (2) 

with Winv the inversion splitting and /i the dipóle mo- 
ment. For ^^NDs the inversion splitting is 1430.3 MHz 
and the dipóle moment is taken to be the same as for 
"NDg, i.e. 1.48 D 0. At even larger field strengths, 
the effect of mixing with higher rotational levels becomes 
more important and higher order terms have to be taken 
into account. 

Using the four electrode geometry shown in Fig. ^ 
it is possible to create a field with dominant dipolar, 
quadrupolar or hexapolar contributions, by applying the 
appropriate voltages. The geometry of the trap is the 
same as was used previously to create an ac electric trap 
|23j |. It consists of two rings with an inner diameter of 
10 mm and two endcaps. The two rings have a width 
of 4.6 mm and are separated by 2.9 mm, whereas the 
two endcaps are separated by 2zq = 9.1 mm. They each 
have a hole with a diameter of 2.0 mm, the hrst to allow 
molecules to enter and the second to cxtract the ions that 
are produced in the detection schcme. 

In Fig.[21the voltages needed for two different trapping 
fields are shown. When a voltage —U2 is applied to the 
two ring electrodes and a voltage U2 is applied to both 
endcaps, such as shown on the left-hand side of the figuře, 
a field is formed with a dominant quadrupole term. Using 
a commercially available finite element program [23| the 
potential of this configuration is simulated. Below the 
trap cross-section the electric field strength as a function 
of position is shown for both the p-direction (solid line) 
and the z-direction (dotted line) . Fitting the correspond- 
ing potential to equation^ we find that $0 ~ —0.43 ř72, 
$2 ~ 1-05 [/2, $4 « 0.35 U2, and $6 « -0.04 U2 and that 
the uneven terms are zero due to symmetry. Beneath the 
electric field strength in Fig. |21 the Stark energy in mK 
(T = Wstark/k^ with k the Boltzmann constant) of a 
^^NDa molecule in this field is shown. The electric field 
in the quadrupole trap increases linearly with distance, 
but, because of the nonlinear dependence of the Stark en- 
ergy on the electric field, the potential is harmonie near 
the center of the trap. Further away from the center the 
potential is linear. The resulting trap is tight (the trap- 
ping frequency is Wz/27r = 1.7 kHz at the center of the 
trap) and has a depth of 0.09 K for i^NOg. 

A hexapole trap results when the voltages and —Uz 
are applied alternately to the four electrodes of the trap 
as shown on the right-hand side of Fig. [3 A fit of the po- 
tential of this field to equation^ylelds <í>i = O, $3 ~ 
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FIG. 2: Two cross-sections of the trap. On the left-hand side, 
the voltages that result in a quadrupole trap are shown, along 
with lineš of constant electric field for electric field strengths 
of every kV/cm from 1 to 7 kV/cm and every 5 kV/cm from 
10 to 25 kV/cm (shown in color online, with increasing field 
strength going from blue to red). On the right-hand side, 
the same is shown for a hexapole trap. Below the two cross- 
sections, the electric field is shown along with the Stark energy 
of a ^^NDa molecule as a function of position in the p- and 
z-direction (solid line and dotted line, respectively). 

and <í>5 « —0.04 U^, with the even terms being zero due 
to symmetry. Again, both the electric field strength along 
the p- and z-direction (solid and dotted line, respectively) 
and the Stark energy of a ^'"^NDs molecule in this field are 
shown below the cross-sections of the trap. The electric 
field strength in this trap is quadratically dependent on 
the distance from the center. However, due to the non- 
linear Stark effect the potential is essentially fiat near the 
center of the trap. Away from the center the potential 
becomes harmonie. This trap is less tight (the trapping 
frequency is a;z/27r = O at the center of the trap and goes 
to íjjzl'2i^ = 750 Hz further away from the center) and 
has a depth of 0.35 K for ^^NDs. 

When applying voltages of -ř/i, -0.6 ř/i, 0.6 ř/i, and Ui 
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FIG. 3: Two cross-sections of the trap. On either side the 
voltages that are applied when a dipóle field is added to 
a hexapole field are shown. Lines of constant electric field 
are shown for electric field strengths of 0.1 kV/cm, every 
0.5 kV/cm from 0.5 to 2 kV/cm, and every 5 kV/cm from 
5 to 25 kV/cm (shown in color online, with increasing field 
strength going from blue to red). On the left-hand side (right- 
hand side) Ui and U3 have the same (opposite) sign. In each 
case both the electric field and the Stark energy of a ^^NDa 
molecule as a function of position are shown in both the p- 
and z-direction (solid line and dotted line, respectively). 



consecutively to the four electrodes of the trap, a dipolar 
field is formed. Fitting the potential of this field to equa- 
tionClresults in $1 « 0.88 [/i, $3 = O and $5 « 0.09 Ui, 
with the even terms being zero due to symmetry. Adding 
a dipóle field to a hexapole field results in either a double- 
well potential or a donut potential, depending on the rel- 
ative sign of the dipóle and hexapole term. In Fig. |3| 
both configurations are shown, along with the electric 
field strength and the Stark energy of a ^^NDa molecule 
in this field along the p- and z-direction (solid and dot- 
ted line, respectively). Both the distance between the 
two minima in the double-well trap and the diameter of 
the ring-shaped minimum in the donut trap are depen- 



dent on the ratio between the dipóle and hexapole term. 
In the upper graph of Fig.^the distance between the two 
wells is shown as a function of the ratio lUi/Usl in the 
case of an ideál hexapole and dipóle field (solid line). The 
dotted line shows the same for the diameter of the donut- 
shaped minimum. The crosses and triangles (double-well 
trap and donut trap, respectively) indicate the simulated 
value of this well separation, using the finite element pro- 
gram with U3 — 5 kV and different values of Ui . 

The height of the barrier between the two wells is 
shown in the lower graph of Fig. 01 given in mK for 
^^NDa. In the case of an ideál hexapole and dipóle field 
(solid line), this height is the same for the double-well and 
donut trap and is solely dependent on the value of \Ui\. 
Note that the barrier height is identical to the amount 
of potential energy gained by a ^^NDs molecule in the 
center of onc of the wells when switching to the single- 
well configuration. This can be seen by realizing that 
in the center of the double-well trap, i.e. at the top of 
the barrier, the electric field strength is determined by 
the dipóle term only, as the hexapole field is zero there. 
In the center of one of the wells the hexapole field is 
equal to the dipóle field, although opposite in direction. 
Therefore, when switching the dipóle field off to produce 
a single- well trap, molecules originally at the center of 
one of the two wells will now have a potential energy 
equal to the barrier height. Values for the barrier height 
were obtained using the finite element program as well 
and are shown as crosses (double-well trap) and trian- 
gles (donut trap) in the lower graph of Fig. ^ From 
these two graphs it is clear that both the well separation 
and the barrier height, and therefore also the coUision 
energy of the molecules upon double- to single-well con- 
version, can easily be changed by varying the voltages on 
the four electrodes. At the maximum value of |ř7i| used 
in the experiment, |ř7i| = 0.55 kV, the barrier height is 
on the order of 5 mK. The separation between the wells 
is 1.8 mm in this case {U3 — 5 kV), whereas the diameter 
of the donut-shapcd minimum is 2.2 mm. 



III. EXPERIMENT 

In Fig. the experimental setup is shown. In the ex- 
periments, the trap is loaded from a Stark-decelerated 
molecular beam. A mixture of 5 % ^^NDs seeded in 
Xenon is supersonically expanded from a pulsed valve 
at a 10 Hz rate. Due to cooling of the valve to — 70°C, 
the gas pulse has a velocity of 280 m/s. In the expan- 
sion, about 60 % of the ammonia molecules are inter- 
nally cooled to the ground statě of para-ammonia, the 
\J,K) = |1,1) inversion doublet in the electronic and 
vibrational ground statě. After passing through a skim- 
mer, those molecules that reside in the low-field seek- 
ing hyperfine levels of the upper component of the inver- 
sion doublet are focused into the Stark decelerator by a 
6 cm long hexapole. Molecules in high-field seeking lev- 
els are attracted by the high electric fields of both the 
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FIG. 4: Well separation as a function of \Ui/U3\ and barrier 
height (for ^'^NDs) as a function of |ř7i| in the ideál double- 
well (solid line) and donut trap (dotted line). In both graphs 
the crosses (double- well trap) and triangles (donut trap) show 
the same quantities obtained from electric field simulations of 
the actual trap geometry. 
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FIG. 5: Experimental setup. A pulsed supersonic beam of 
^^NDa molecules seeded in Xenon exits a cooled valve with a 
velocity of 280 m/s. Ammonia molecules in the low-field seek- 
ing hyperfine levels of the | J, K) = |1, 1) statě are focused into 
the decelerator by a hexapole, decelerated, and then transver- 
sally and longitudinally focused into the trap by a hexapole 
and buncher, respectively. The molecules are detected in a 
TOF-mass spectrometer. 



hexapoles and the decelerator and will be lost from the 
beam. Molecules that reside in the remainder of the hy- 
perfine levels of the inversion doublet, that only have a 
higher order Stark efFect, will hardly be effected by the 
electric fields. The decelerator consists of 95 pairs of 
parallel 3 mm-diameter rods, with a closest distance of 
2 mm between two rods and a closest distance of 2.5 mm 



between pairs of rods. Each consecutive pair of rods is 
rotated over 90 ° and within each pair, a positive voltage 
of 10 kV is applied to one electrode, whereas the same 
negative voltage is applied to the other. A more extensive 
description of this part of the setup, along with the oper- 
ati on p rinciple of the decelerator can be found elsewhere 
[23. Using a phase angle of 57.5° a subset of low- 
field seeking molecules with an initial velocity of around 
280 m/s is decelerated to about 15 m/s at the exit of 
the decelerator. After 5 mm the molecules are trans ver- 
sally focused into the trap by a second (12.5 mm long) 
hexapole. Subsequently, they are longitudinally focused 
into the trap using a buncher that is placed another 5 mm 
downstream of the hexapole. The buncher has a similar 
geometry as the trap. Before the ammonia molecules en- 
ter the trap, voltages are applied to the trap such that a 
last electric field slope is formed. In gaining Stark energy 
when entering the trap, the low-field seeking molecules 
lose their last kinetic energy and reach a stand-still in 
the middle of the trap. At this moment the voltages on 
the four electrodes are switched to generate a particular 
trapping potential, thereby confining the molecules in ei- 
ther a quadrupole, hexapole, double- well or donut poten- 
tial. After a certain trapping time, the trap is turned off 
and molecules in the upper component of the inversion 
doublet are ionized in a state-selective (2+l)-Resonance 
Enhanced Multi Photon lonization (REMPI) scheme, us- 
ing a pulsed UV-lascr around 317 nm 26]. The rcsulting 
ions are subsequently detected in a time-of-flight (TOF) 
mass-spectrometer. Using a 75 cm lens the ionization 
laser is focused into the trap passing through the 2.9 mm 
gap between the two ring electrodes. The laser focus is 
estimated to be about 200 /im in diameter. Applying 
small bias voltages of 200 V, 150 V, -150 V and -200 V 
to the four electrodes of the trap results in an electric 
field that accelerates the ions towards the detector. The 
ions are detected over a length of 2 mm along the laser- 
beam, limited by the size of the hole in the endcap. Us- 
ing a Computer controUed translation stage, the focus of 
the laserbeam can be moved along the beam-axis, i.e. 
along the z-direction, over about 2.5 mm. In this way, 
the spatial distribution of the trapped molecules can be 
measured. 



IV. RESULTS AND DISCUSSION 

In Fig. El the density of trapped ^^NDs molecules is 
shown as a function of position in the z-direction, both 
for hexapole and for quadrupole trapping voltages. The 
measurements are taken after the molecules have been 
trapped for about 72 ms. We found that in order to 
obtain a single-well profile for the molecules confincd in 
the hexapole trap a small dipóle term had to be added. 
We attribute this to deviations of the distances between 
the electrodes from the ideál geometry. This compen- 
sating dipóle term is added in all subsequent measure- 
ments. As the quadrupole trap is much fighter than the 
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FIG. 6: Density of trapped ^^NDs molecules as a function of 
position in the z-direction for both the quadrupole trap (upper 
trace) and the hexapole trap (lower trace), after trapping for 
about 72 ms. 

hexapole trap (see Fig. [SJl, a signál difference arises be- 
tween the two traces. The peak density of molecules in 
the hexapole trap is about lO^cm^'^, whereas it is about 
ten times higher in the quadrupole trap. The loading 
process is identical for both traps and is matched to the 
hexapole trap, resulting in a less than perfect match of 
the decelerated packet to the quadrupole trap. This gives 
rise to the peculiar shape of the measured distribution in 
the quadrupole trap, which seems to consist of two pack- 
ets with different temperatures. 

In Fig. [7| the density of ^^NDa molecules in the 
hexapole trap is shown as a function of position in the 
z-direction. By adding a dipóle term to the hexapole 
term, the shape of the trap and thereby the profile of the 
trapped molecules can be changed significantly. Each 
trace in Fig. corresponds to a different value of the 
added dipóle term. When no dipóle is added, the trap 
consists of only a single well. A single peak results, as 
trace (d) shows. Adding a negative dipóle term trans- 
forms the single well to a double well, with a distance and 
barrier between the two wells that depend on the value 
of the dipóle term. The upper three traces of Fig.[7|show 
two peaks that move closer together as the dipóle term 
is increased from -0.55 kV (trace (a)) to -0.3 kV (trace 
(b)) to -0.15 kV (trace (c)). When a positive dipóle term 
is added, the shape of the trap becomes that of a donut. 
As only molecules in the relatively small laser focus ncar 
the center of the trap are ionized in the detection scheme, 
this transformation to a donut-shaped trap leads to a re- 
duction in signál only, as is shown by trace (e) and (f) 
{Ui = 0.25 kV and Ui = 0.55 kV, respectively). 

At the densities that we have in the trap, no thermal- 
izing coUisions occur, and formally no temperature can 
be assigned. To nevertheless characterize the trapped 
sample in terms of a temperature, we fitted the density 
profile to the thermal density distribution p^ : 

n{p, z) = „oe-^«""-'=('''^)/'=^^, (3) 

with no the density at the trap minimum and 
Wstark {p, z) the potential energy as a function of posi- 
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FIG. 7: Density of ^^NDa molecules confined in a hexapole 
trap as a function of position in the z-direction. A negative 
dipóle term of -0.55 kV, -0.3 kV and -0.15 k V is added to the 
hexapole term (ř/s = 5 kV) in the upper three traces (a-c), 
resulting in a double-well potential. Trace (d) shows a packet 
in a purely single well. A positive dipóle term of 0.25 kV and 
0.55 kV is added for the lower two traces (e-f), thereby creat- 
ing a donut potential. On the right-hand side, trace (a) and 
(d) are shown once more, along with simulations of the den- 
sity distribution of the packet of molecules for temperatures 
of 5, 10 and 20 mK. 



tion. On the right-hand side of Fig. [7| the measurements 
with a dipóle term of Ui = —0.55 kV and Ui = ú kV 
are shown once more, along with the thermal density 
distribution as obtained from equation ^ for tempera- 
tures of 5, 10 and 20 mK. It can be seen that in either 
case the simulated distribution with a temperature of 
10 mK matches the measurement best. The same sim- 
ulation has been performed for the spatial distribution 
of molecules confined in a quadrupole trap. Due to the 
seemingly double profile, it is more difhcult to describe 
this distribution with a single temperature. A temper- 
ature of 75 mK seems to best describe the broad base, 
whereas the sharp center peak agrees best with a tem- 
perature of 5 mK. An alternativě method to deduce a 
temperature is to measure the width of the packet of 
trapped molecules as a function of time after the trap 
has been turned off. Measurements of this type are 
shown in Fig. |H1 On the left-hand side two profiles of 
^^NDa molecules are seen, one after 20 /is and one after 
420 /KS after switching off the hexapole trap (upper and 
lower curve, respectively). In either case the molecules 
have been confined for about 72 ms. Due to the veloc- 
ity spread of the molecules the packet spreads out af- 
ter the trap has been turned off, resulting in less signál 
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FIG. 8: Determination of the temperature of a packet of 
^^NDa molecules trapped in the hexapole trap. On the left, 
the profile of a paclíet is shown 20 /xs and 420 /is after the 
trap is turned off (upper and lower trace, respectively). From 
a fit to a Gaussian lineshape the FWHM of the distribution 
is obtained, as shown as a function of free-fiight time (solid 
line with crosses) on the right-hand side of the figuře. 
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FIG. 9: Density of ^^NDa molecules in the quadrupole (upper 
trace) and the hexapole (lower trace) trap as a function of 
time that the trap is on. 

and an increased width. To obtain this width a Gaus- 
sian lineshape is fitted to each profile. The acquired 
fuU-width-half-maximum (FWHM) is shown as a func- 
tion of free-flight time on the right-hand side of the fig- 
uře (solid line with crosses). To deduce a temperature, 
this curve is fitted to \J (Az)^ -|- (Aw^ • í)^ (dotted line 
with triangles), resulting in a velocity spread of 2.3 m/s 
and a position spread of 1.8 mm (both FWHM). Using 
Au^ = V81n2 • yJkT/m El], this gives a temperature of 
about 2.5 mK. 

The different temperatures of 10 mK and 2.5 mK ob- 
tained by the two different methods actually correspond 
to a difference of only a few tenths of a mm in the width 
of the initial packet. Such a small difference can be ex- 
plained by possiblc misalignmcnts of the electrodes of the 
trap, causing uncertainties in the actual electric field, and 
therefore in Wstark{p, z). Moreover, extracting the width 
of the packets from the measurements is quite sensitive 
to how the baseline is taken into account in the fit. All 



together, this analysis only serveš to indicate that the 
temperature of the trapped sample is in the mK range. 

In Fig.Elthe density of ammonia molecules is shown in 
both the quadrupole and the hexapole trap as a func- 
tion of time that the trap is on. The molecules are 
seen to leavc the trap with a 1/e time of 0.16 s for the 
hexapole trap and 0.17 s for the quadrupole trap, the 
same within error bars. The trapping lifetime is mainly 
limited by background coUisions (background pressure 
4- 10~* mbar). Optical pumping by black-body radiation 
cannot be completely neglected, however. In the room- 
temperature trap, pumping of the ^^NDs molecules to 
the vibrationally excited V2 level (umbrella mode) around 
750 cm^^ is expected to occur on the time-scale of a sec- 
ond. 



V. CONCLUSIONS 

We have presented a new electrostatic trapping geom- 
etry, that illustrates the versatility of electric fields in 
forming tailored trapping potentials. By applying differ- 
ent voltages to a set of four electrodes, different trapping 
fields are created, between which rapid switching on a 
timescale of 100 ns is possible. Spatial distributions and 
lifetimes of a cold packet of ^^NDs molecules loaded into 
both a hexapole and a quadrupole trap are shown. Ad- 
ditionally, we studied a hexapole trap that transforms 
into a double-well or donut trap when different dipóle 
terms are added. The barrier height and well separation 
of either one of these two traps depend on the strength 
of the dipolar and hexapolar electric field and are easily 
changed. 

Rapid switching between, for instance, the double-well 
trap and a pure hexapole trap offers good prospects for 
measuring collision cross-sections as a function of col- 
lision energy at low temperatures. The amount of en- 
ergy gained by molecules trapped in one of the wells 
after switching to the single-well depends solely on the 
strength of the electric dipóle field, which can be easily 
varied. 
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